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Abstract 

We propose a dynamical model for the unpolarized generalized parton 
distribution function (GPD) based on impact parameter representation and 
discuss its x and t dependencies at large and small values of these variables. 
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Recent developments in the nucleon structure studies are closely related to the 
generalized parton distributions (GPDs) H q (x,£,t) and others. These functions 
are the matrix elements of the quark field operators between nucleon states of the 
different momenta with x = l/2(xj + Xi) and £ = 1/2 (x/ — Xj)Q. This topic 
attracted much attention in the last decade. The results are collected in the re- 
view papers flU |2l [3l |H [51 [6l 13 and many references to the original papers can be 
found therein. The important point is that the Fourier transform of the function 
H q (x, £, t) can be related to a weighted parton number density corresponding to 
the quark q flavor in the transverse plane when interpolated to the point £ = 
where the variable £ is the longitudinal momentum transfer ( skewedness pa- 
rameter). The impact parameter-dependent GPD allows one to obtain information 
on the parton distribution in the transverse plane over the nucleon area simultane- 
ously with the parton distribution over fraction x of the total nucleon momentum. 
This result would pave the roads from the experimental data related to the ampli- 
tudes of the processes y*p — > yp (deeply virtual Compton scattering (DVCS)) and 
hard meson exclusive production y*P ~~ * Mp to the studies of hadron structure in 
the two-dimensional transverse space. For the first time position space aspects of 
parton model were studied in the paper flSJ. 

A common procedure of obtaining GPD dependence is to guess the form of 
H q (x,£,t) proceeding from the model assumptions, namely use of the general 
constraints based on form-factors, first principles, experience gained from high- 
energy scattering and experimental studies of DVCS [fT0l[TTl[T2l[T^[T4l[T5l[T6ll . 
However, since impact parameter representation provides a useful framework for 
implementing geometrical aspects as well as transparent account of the unitar- 
ity constraints, it can be used as a complementary way for the GPD's modelling. 
In particular, unitarity of the quark-hadron scattering amplitude at low values of 
x can be implemented in that way [1T8 1 and forms of gluon distributions were 
elaborated on the base of this representation in ffT9ll20l . Impact parameter repre- 
sentation is also very helpful in establishing the relation of the GPD dependence 
with color transparency phenomena [|2Tll22ll23l . In this note we continue to pro- 
ceed that way and construct a model for the impact parameter dependent GPD to 
provide further insight into the spatial hadron properties. 

Prior to the model description we would like to provide several known gen- 
eral relations for GPDs. We consider generalized parton distribution function 
H q (x, £, t), i.e. the case of unpolarized hadron and parton. In the forward limit 
this function reproduces standard, integrated over impact parameter, parton distri- 
bution function, H q (x, 0, 0) = q(x). If one uses the combination 

H q (x, t) = H q (x, £ = 0, t) + H q (-x, £ = 0,t), 

'Note, that Xi and Xf are the light-cone momentum fractions of the quark in the initial and 
final state, respectively. 
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then the nucleon Dirac form factor F 1 (t) is the sum F 1 (t) = J2 q e qF?(t) where 

Fi (t) is represented as the integral Ff(t) = H q (x, t). Relation of the impact 
parameter dependent parton distribution with the corresponding GPD was estab- 
lished in [24] and has the form 



l r°° . . . 

q(x } b) = — s/^tdV^iH q (x,t)J (bV^i), (1) 
27r Jo 

where b = |b| is the distance from the active quark position to the hadron center 
of mass in the transverse plane, which is determined by the relation 



0, (2) 



and Xi stands for the longitudinal momentum fraction of the parton i in the infinite 
momentum frame. It should be noted that the function q(x, b), in contrast to the 
function q(x), is dimensional, it has dimension of a squared mass and is inter- 
preted as a parton number density referred to the transverse area of a nucleon. It 
can be written therefore in the form q(x, b) = fi 2 q(x, b), where fi is a some mass 
scale related, as it will be clarified later, to the finite transverse size of a parton 
and q(x, b) is a dimensionless parton number density distribution, i.e. it can be in- 
terpreted as a probability to find quark of flavor q with the the momentum fraction 
x of the total nucleon momentum at the distance b from the hadron's center. 

Another important parameter is the distance from the active quark to the center 
of mass of other constituents (spectators) 

rx = (3) 
1 — X 

Its mean value can also be associated with the hadron size. Requirement that the 
value (r 2 ± ) x remains finite at x — > 1 leads to t-dependence of GPD in the form 
of the product (1 — x) n t, n > 2 in the limit x — > 1 [|22l . Another requirement 
is that the values of (h 2 ) x and, of course, (r\) x , should not increase faster than 
In 2 - at x — > 0. It has grounds in the short-range nature of strong interactions and 
unitarity. 

Keeping in mind the above preconditions, we propose the following form for 
the function q(x, b): 

q{x,b) = (4) 
1 + q {x,b) 

Such form for q(x, b) can be obtained with account for the finite size of the parton 
in the transverse area ( determined by the inverse mass 1//2) in a way similar to 
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the one used in lfT71 for the case of the gas of extended hadronsH. The form of Eq. 
H]can also be related to the form of the quark-hadron amplitude as it was discussed 
in flU. 

The function qo(x, b) is a dimensionless distribution function for the point-like 
partons. It should be noted that q(x, b) ~ qo(x, b) in the region where qo(x, b) <C 
1. It takes place, e.g. at large values of b, where the finite size of partons does not 
play a role. 

Construction of the function q(x, b) in the form of Eq. |4] implies that the 
weighted parton number density referred to the transverse area of the nucleon 
cannot exceed the magnitude of inverse parton transverse area. This form and the 
model for the function %(x, b) described below lead to the saturation of the parton 
distribution function q(x, b) at low values of x and b. The mass scale pt, determines 
transverse parton size and can be related to the constituent quark radius. We note 
also, that the following upper bound takes place: 

q(x,b) < 1/Sg, 

where S q is a transverse area of an extended parton. Indeed, there are experimental 
indications on the presence of the extended objects inside the proton found in the 
DIS data 1251. 

To construct the function q~o(x, b), a dimensionless distribution of the point- 
like partons, we consider the picture of hadron structure described in [27 J. In 
this approach hadrons contain massive quarks interacting with self-consistent pion 
field and the pions themselves are not elementary. They consist of strongly bounded 
quark- antiquark pairs. This picture is a model-independent consequence of the 
spontaneous symmetry breaking. Thus, the active quark moving in the field of the 
spectators interacts with them by a pion exchange. This is strong interaction and 
it cannot be reduced therefore to one-pion exchange, and spectators should not be 
considered as the remaining valence quarks only. Being essentially relativistic, 
the bound state problem should be addressed on the multiparticle basis and Fock 
components with many partons should be taken into account. Thus, the dynam- 
ics of the quark interactions and probability to find quark at the distance b from 
the hadron center depend, in fact, on the distances between active and spectator 
quarks. After averaging over positions of the spectator quarks, the r_i_ = \r±\ will 
be the only distance the function %(x, b) should depend on. The described mech- 
anism implies that the lifetimes of the active and spectator quarks are comparable, 
i.e. the process is coherent. Such a coherence can occur even in the case when 
an active quark is hard and has large value of x and spectator quarks are soft and 

2 Indeed, the number density of point-like partons no can be written in the form no = 
(s h -s N) > where N is the total number or partons (located in the hadron transverse area Sh) 
with transverse size S q . 
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have small values of x. The interaction of the active quark with spectator ones 
should be considered as a final-state process as it was discussed in [28]. 
For the input distribution function qo(x, b) we use a factorized form 

With account for quark interactions with the pion field a simple exponential de- 
pendence was chosen for the function p(r±), 

p(r ± ) ~ exp(-fxrj_), (5) 

where mass scale p is order of the pion mass (but should not be exactly equal to 
it). Of course, this is oversimplified form because of unspecified pre-exponential 
factors dependent on the variable r± can be presented also. For the purposes of 
qualitative consideration we neglect these complications. The function xi x ) can 
be connected somehow to a number of the spectator quarks interacting with active 
quark (having the fraction of nucleon momentum equal to x). 

For the behavior of x{ x ) m tne limiting cases x — > and x — > 1, we use the 
standard dependencies borrowed from the deep-inelastic processes studies of the 
parton distribution functions: 

x(x) ~ x~ x and x( x ) ~ (1 — x) (6) 

respectively. Note, that parameter A > 0. Thus, the adopted expression for the 
function qo(x, b) is the following 

Qo(x, b) = x(x) exp(--^-). (7) 

1 — x 

The function x( x ) can ' m principle, be obtained from the parton distribution 
q(x) with the use of Eq. ©, i.e. due to relation 

g(x)~-(l-x) 2 Li 2 [- X (x)], (8) 

where the function Li 2 (z) is dilogarithrrjl Qualitative behavior of the dimension- 
less parton distribution q(x, b) with the function x{ x ) having limiting dependen- 
cies in the form of Eq. © is illustrated by the plot in Fig. 1. Evidently, the 
function q(x, b) is saturated in the region of small x and b. 
The average value (h 2 ) x is determined by the relation 

= JiTj^gfo b ) bdb (9) 

fo° <l{ x y b)bdb 



3 The dilogarithm is defined by the series hi2(z) — X^fcLi w anc ^ at l z l — 1 tn i s function 
defined through an analytic continuation. 
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Figure 1 : Dimensionless generalized parton distribution function q(x, b) (color 
online). 

and 

It is clear that q(x, b) at large values of x has the same dependence as the function 
qo(x, b) and it immediately provides at x — > 1 the finite value for (r 2 )^, i.e. 

{r 2 ± ) x ~ 1/yU 2 , at X -»• 1 

while 

(b% ~ (1 - x) 2 //i 2 . 

At a; — >• the expression for the function q(x, b) in the form of Eq. © with 
limiting dependence of Eq. © provides 

(b%-(ri)^|ln 2 (i). 

Thus, the model leads to reasonable results for the hadronic sizes at small 
and large values of x. At large values of x the hadron size is ^-independent. 
It has finite value and consistent with quark confinement. At small values of x 
hadronic size increases in a way consistent with unitarity and short-range nature 
of strong interactions. Such behavior is a result of transform of the x -independent 
6-dependence of qo(x,b) (at small values of x) into the function q(x,b) which 
results in average squared impact parameter 

(b 2 ) x = 4^-lnH(x,t)\ t=0 , 
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rising logarithmically with x — » 0, i.e. a/ (b 2 ) x ~ m (i)- ^ ne can concm de 
therefore that hadrons participating in hard interaction have finite radius (and this 
does not imply color transparency phenomenon), while soft hadron interaction 
radius is increasing with energy like In s (since s ~ 1/x). 

We consider now the x and t-dependencies of the function H(x,t) at large 
values of x: 

"' (M) ~ [l-(l-s)WF 2 ' (10) 
This form is the Fourier-Bessel transform of the function q(s, b). From Eq. (fTOl) it 
is evident that at fixed values of —t the function H q (x, t) does not depend on this 
variable in the limit x — > 1. If we take the fixed values of x and large values of 
—t the function Tt q (x,t) seems to decrease as (— t)~ 3 / 2 . But it is not legitimate 
to consider the above two limits independently, separating x and t dependencies, 
because those limits are not consistent simultaneously with the requirement for 
the parton process to be coherent, which is a precondition for the validity of the 
discussed mechanism. Instead, we should impose additional relations which guar- 
antee coherence at large values of —t, e.g. 

(1 — x)t = const. (11) 

or 



(1 — x)v— t = const. (12) 

Eq. (fTTj) is known as the Berger-Brodsky limit ll29ll and it provides the condition 
for the process of parton interactions remain coherent as well as Eq. (PT2]) does (cf. 
[|28l ). Thus, when we consider asymptotical dependence on —t we should always 
take x — > 1, otherwise coherence would be destroyed. Asymptotical behavior of 
form-factor, which is an integral 



F?(t) = [ dxH q {x,t), 
Jo 



has the following power-like form : 



N+3 
2 



In the region of small values of x the function TC q (x,t) has the following 

w 



dependencej at small —t: 

H g (x,t) ~ In 2 ^— ^ exp 



A 2 , 2 (\ 
:t In 2 - 



V x 



(13) 



4 It results from behavior of the Bessel function Ji (x) and x( x ) at small values of x. 
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In the framework of the adopted picture for the hadron structure this low-x de- 
pendence should be associated with the sea quarks. Eq. (fT3l) is in a qualitative 
agreement with the experimentally observed exponential t-dependence of the dif- 
ferential cross-section of the DVCS process lf30l . 

Thus, we can conclude that the GPD is no constrained by the form factor 
alone, which is sensitive to the difference of quark-antiquark number densities 
only. Instead, in addition to the amplitudes of the DVCS processes, it would be 
promising to relate GPD behavior to the amplitudes of another exclusive process 
such as elastic hadron-hadron scattering and use experimentally measured differ- 
ential cross-section of elastic scattering to constrain the t-dependence of GPDs. 
Large angle elastic scattering is sensitive to the collisions with small impact pa- 
rameter and provides information on the hadron structure in the central region in 
the position space. 

Another possibility is related to the multiparticle production. Proceeding along 
the lines described in OTA , it is possible to treat the amplitude of such process 

F(Bi - B 2 ,bi,x 1 ....b n ,x n ,b' 1 ,x , 1 ....b' n ,,x' n ,) 

as the wave function of the production of the two clusters with n and n' parti- 
cles with centers located in the position space at Bi = Y^h=i x i^i and B 2 = 
Y^i=i x 'Mi- Using idea of quark-hadron duality, one can try then to extract in- 
formation on the GPDs from the observables measured in the hadron production 
processes. This possibility will be pursued elsewhere. Here we note that the 
hard inclusive processes are complementary to hard elastic scattering since they 
are sensitive not only to the collisions with small impact parameters, but to the 
peripheral collisions also [32]. 

To conclude, we note that the proposed model for the impact parameter-dependent 
GPD allows one to consider simultaneously the regions of large and small values 
of x and —t. The parton distribution q(x, b) is saturated at small values of x and 
b. It provides finite value for the hadron size determined by the inverse pion mass 
at large values of x and logarithmically growing hadron size at small values of x. 
This growth is due to saturation of the weighted probability distribution q(x, b) in 
the impact parameter representation. 

We are grateful to A. Kisselev, V. Petrov and B. Pire for the interesting discus- 
sions and useful remarks on the list of references. 
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